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Abstract 

! We find explicit expressions for two first finite size corrections to the distribution of Bethe roots, the 

asymptotics of energy and high conserved charges in the st(2) quantum Heisenberg spin chain of length J 
in the thermodynamical limit J —> oo for low lying states with energies E ~ 1/J. This limit was recently 
studied in the context of integrability in perturbative Af = 4 super- Yang-Mills theory. We applied the double 
Qh! scaling technique to Baxter equation, similarly to the one used for large random matrices near the edge of the 

eigenvalue distribution. The positions of Bethe roots are described near the edge by zeros of Airy function. 
Our method can be generalized to any order in 1/J. It should also work for other quantum integrable models. 
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1 Introduction 



We study in this paper the integrable periodic Hciscnbcrg XXX S chain of noncompact quantum spins 
transforming under the representation s = —1/2 of sl(2), in the so called thcrmodynamical limit of large J, 
where J is the length of the chain, in the ferromagnetic regime of low energies 

The problem is known to be solvable by the Bethe ansatz approach (see for example [I]) and the energy 
of a state of S magnons in dimcnsionlcss units is given by a simple formula 



where the Bethe roots Uj, j = 1, 2, • • • , S, parametrizing the momenta of magnons, are solutions of a system 
of polynomial Bethe ansatz equations (BAE) 



It can be proven that for this model the roots are always real. 

In the thermodynamical limit we will also consider a large number of magnons S ~ J. It is clear then from 
eq. (TTJ) that in order to focus on the low lying states with energies E ~ 1/ J we should take the characteristic 
Bethe roots of the order Uj ~ J. It means that the chain is very long and the spins are rarely changing 
along it. The typical length of spin- waves (magnons) is of the order of the length J. Our goal is to study 
the limiting J — > oo distributions of Bethe roots and the finite volume 1/ J corrections to these distributions, 
to the energy and higher conserved charges. In the main order this thermodynamical limit for the compact 
Heisenberg XXX 1 / 2 chain of su(2) spins was already considered in [5], and later in [3] in relation to the 
integrable dilatation hamiltonian in planar perturbative superconformal AT = 4 super- Yang-Mills (SYM) 
theory. Its description and the general solution in terms of algebraic curves was proposed in [3] for the 
su(2) cas<3 and in [51 [7] for the sl(2) chain. We will concentrate in the current paper on this last case. 
Generalization to the su(2) case is straightforward. 

The study of 1/ J corrections in these systems was started recently in the papers [HI H] for the simplest 
single support, or one cut distribution, whereas a similar quantum h correction to the classical KdV solitons 
was already found earlier in the general multi-cut case in [10j . 

The main results of our paper are: 

1. The explicit formulas for the 1/J and 1/J 2 corrections to the general multi-cut distribution of Bethe 
roots and to the corresponding energy of a Bethe state in terms of the underlying algebraic curve. 

2. The universal description of the distribution of Bethe roots in the vicinity of an edge of a support in 
terms of zeroes of the Airy function, similar to the double scaling limit in the matrix models. 

3. Asymptotics of conserved local charges Q n (S, J) in the large n limit. 

4. Asymptotics of conserved global (non-local) charges of high order. 

Unlike the papers [S] |H] using the method of singular integral equation corrected by so called anomaly 
terrrQ we will use here the exact Baxter equation written directly for the analytical function - the resolvent 
of the root distribution (similar approach was used in [B]). This equation is valid before any limit. 

These results might be interesting for different kinds of specialists. 

First, for those who are studying large N random matrix or random partition ensembles. In particular, 
the distributions of eigenvalues in matrix models are described by similar (although not the same) algebraic 

It is different from a more traditional regime E ~ J widely studied since many years, especially in the condensed matter 
literature. 

2 Following a similar approach of [5] to a somewhat different limit of large spin 

3 This phenomenon of anomaly, or the contribution of close eigenvalues in the thermodynamical limit of BAE was first 
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observed in [IT] 
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curves. The 1/J corrections remind the l/N corrections, and the Airy edge distribution observed in our 
work is known to describe the edge behaviour in the double scaled matrix integrals as well [12]. No doubt 
that by choosing particular solutions of cq.Q or modifying its l.h.s. (by considering inhomogencous chains) 
we can also find here various multi-critical phenomena similar to those found in the matrix models in the 
context of quantum 2D gravity applications |13j . 

Second, it might be interesting to those who work on various aspects of AdS/CFT correspondence in 
supersymmetric string and gauge theories and in particular on the integrability in Af = 4 SYM theory (see 
[55] and references therein). They might shed some light on quantum corrections to the classical limit of the 
AdS dual of SYM, the superstring on the AdS§ x S§ background, now known only for particular classical 
solitonic solutions [TS] of rotating strings. AdS/CFT correspondence should manifest itself as the coincidence 
of such corrections in these so different intcgrable systems. Their similarity, and even the coincidence in a 
certain regime, was already observed on particular string and chain solutions, having only one support for 
the Bethe roots distribution [E1IH1IMEDIIBIISESII2- 1/J corrections were first studied for BMN states 
in [3D], where the integrable spin chain for Af = 4 SYM was first proposed, and then in [51]. The Airy 
edge behavior also seems to be universal enough to manifest itself on both sides of AdS/CFT duality. This 
behaviour observed in the present paper for spin chains describing the spectrum of anomalous dimensions 
in the perturbative SYM in thermodynamical limit, is natural to expect also in the quantized string theory 
where the classical limit, the analogue of the thermodynamical limit, is described by similar algebraic curves 
obtained by the finite gap method @. The universal edge behaviour might be completely driven by the 
similarity of integrable structures of these two seemingly different systems. 

The paper is organized as follows. In section 2 we summarize the explicit formulas for the Hamiltonian 
and the local and non-local conserved charges of the model. In section 3 we review the general solution of 
the model in thermodynamical limit when J — > oo and E ~ 1/J. In section 4 we calculate the first 1/J 
correction to the general multi-cut solution from Baxter equation. In section 5 we apply the method of the 
double scaling limit to the Baxter equation and find the near branch point behavior. In section 6 we combine 
the results of section 4 and 5 to find the explicit formula for the second, 1/J 2 correction. In section 7 wc 
give the formulas for 1/J corrections to the energy and compute the asymptotics of high conserved local and 
non-local charges in all orders in 1/J. Conclusions are devoted to unsolved problems and perspectives, as 
well as to the discussion of parallels with the matrix models. In appendix A we write explicit formulas for 
1/J expansion of B AE and in Appendix B we express 1/J 2 corrections to the energy for one-cut solution 
through certain double sums in mode numbers, generalizing the results of [§]. 



2 Hamiltonian, Transfer-matrix and Higher Charges of s[(2) chain 

The hamiltonian of interaction of the neighboring spins si, s;+i can be written in an explicit way [14j 

,/ 

H-l/2 = X] H -l/2 ( 3 ) 

1=1 

with the Hamiltonian density 

H l J l 1 + / 1 2 \k,m-k) = jr (8 k ^{h{k) + h{m-k))-^j-\\k',m-k'), (4) 



k'=0 



where \ki, . . . , fc;, fcz+i, . . . , kj) is a state vector labeled by J integers kj (s = —1/2 spin components) and 
h(k) = yij— i j axe harmonic numbers. 

The total momentum P(u) 

~ u 3 + i/2 1 ' 



See 1111 1201 for this approach on both sides of duality 
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satisfies the (quasi-)periodicity condition following directly from eq.(H]) 

s 

Ptot = p K) = 27rk / J > keZ - ( 6 ) 

In some applications, such as the spectrum of anomalous dimensions of operators in TV = 4 SYM theory, 
we select only purely periodic Bcthc states 

P tot = 27rm, m € Z. (7) 

We can also study other physically interesting quantities of this model, such as the local conserved charges 
Q r . They are defined as follows 

f(v) =cxp^f]Q r ^- 1 ^ , (8) 

where the quantum transfer matrix T(v) = T(v; 0, 0, • • • , 0) is a particular case of the inhomogeneous transfer 
matrix 

T(v,V!,--- ,vj) =Tr Ro,i(v - i>i) • • --Ro,j(v - vj) (9) 



and Rqj is the universal st(2) R-matrix defined as [5T] 



v — ik 



R s(v) =y, r m-p { q % = n (io) 

j=0 fe=l 

with Vqi being the operator projecting the direct product of two neighboring spins sq = si = — 1/2 to the 
representation j . Recall that 

T(v,V!,- ■ ■ ,Vj),T(v';vi,- ■ ■ ,vj) =0 (11) 

for any pair v,v f , due to Yang-Baxter equations on the i?-matrix. 

The direct calculation shows that Ptot = —Qi is the operator of the momentum, such that and H_i/ 2 = 
(§2 is the hamiltonian cq.©, etc. Those charges are local, in the sence that the charge density of Qk contains 
< k consecutive spins. 

There is a more efficient way to generate higher charges than using eq.©. We can use the recurrence 
relation (see [T4"]h 

rQ r+1 = [B,Q r ], (12) 

where B = i^f—il H l -i/2 is the boost operator. The last formula is true up to some boundary terms 
destroying the periodicity, which should be dropped. 

Due to the integrability manifestly expressed by eq. (|ll|) all these charges commute and their eigenvalues 
on a Bcthe state characterized by a set of Bethe roots satisfying eq. ((2|) (enforcing the periodicity of the chain 
or the quasi-periodicity of the Bethe state) are given by [22] 

Qr = ^ ^ ( + */ 2 ) r ~ 1 ~ - l / 2 ^ ) ' (13) 

We can also study the non-local charges. They can be defined in many different ways. The definition 
could be similar to eq.®, but the expansion of T(u) goes around u = oo. However, the most natural charges 
are defined through the resolvent of Bcthe roots: 

S oo 

G(x)=J2—^— = J2d n x- n . (14) 



5 The operators of the type Tr (V fcl Z ■ ■ • V^-'Z) in SYM, where V = d + A is a covariant derivative in a null direction and 
Z is a complex scalar, represent the state vectors \ki, . . . , fc;, fei+i, . • • , kj) and the dilatation hamiltonian is given at one loop 
by the XXX_i/2 hamiltonian. 
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The j-th charge is the j-ih symmetric polynomial of Bethe roots 

i & _i 00 

fe=l k=0 

We will later estimate the behavior of d k ^ n at n — > 00 and high orders of 1/J expansion. This asymptotics 
will be universal and a few leading terms of it will be the same for various definitions of non-local charges. 



3 1/J expansion of B AE 



Let us start from reviewing the "old" method of solving eq.([2]) in the thermodynamical limit J — » 00, 
Uk ~ J ~ S", before sticking with the most efficient one using the Baxter equation. 

As we mentioned in the introduction the eq.© has only real solutions, i.e. all the roots lie on the real 
axis. We label the roots so that ttj+i > Uj. Suppose there exists a smooth function X(x) parametrizing the 
Bethe roots 

u k = JX(k/J), g(X(x)) ee -1- ~ . (16) 

X'(x) Uk+i - u k 

For large S the function g{x) has a meaning of density of Bethe roots. As follows from definition |[T6]) its 
normalization is 

dxg(x) = a (17) 



with a = S/J. Taking log of both parts of eq. @ we havc@ 



u - — i/2 S ' 

2mm j + J log — — - = — 2i arctan(uj — Uk), (18) 

J 1 fe=i 

where rrij are strictly ordered integers m J+1 > rrij. In terms of the logarithm we can write 

u i — i/2 v-V, Uj — Uk + i 
2mnj + J log 3 ! = V log ^ 5 : , (19) 

where rij = rrij — j + arc non-decreesing integers 0. 

Now we have instead of "fermionic" m^'s the "bosonic" rij's which are simply ordered %+i > rij and 
different Bethe roots can have the same magnon numbers rij = rij+i = • • ■ . In the thermodynamical limit 
we can rewrite ea.(fT9f assuming k to be far from the edges, as follows (see also Appendix A) 

E' fu -u k +i\_ „v-v 1 2 V^' 1 2 Y^' 1 2 Y^' 1 

og \u 3 - Uk -i) = ^^~u~ k +T i^ (u 3 -u k Y~l^ ( Uj - Uk) 5 + 7Z-. ( Uj - Uk y (20) 




n q' [coih(ir g)]e 



o 2 / // 

— Z7T p g 



1 



sinh(7T£)) 



7rg (3 '[coth(7rg)]4 



were we introduce the notation defined by [f(g)] n = f(g) — J27=o /^(0)t f° r the functions regular at zero. 
For singular functions the Taylor series should be substituted by the Laurent series so that [f(g)] n is zero 
for gi = and has first n — 1 zero derivatives at this point. The terms in the first line represent the naive 
expansion of the l.h.s. in l/(uj — u k ). It works well for the terms in the sum with Uj S> u k . The terms in 



6 Note that i log 2±» = arctan(x) — ^-sign(a:) for standard definition of the log i.e. log(x) = — log(a: — iO) for x < 0. 
7 This fact is obvious when S <^ J, for J — > 00, when we can neglect the r.h.s. of ea. l|19| l. We believe that this is true also 
in general, although this fact is irrelevant for the rest of the paper 
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POO 









• -1- 

•x=0 






Figure 1: Hyperelliptic Riemann surface 



the second line describe the anomalous contribution at Uj 
can expand 

j - k 



Uk, for close roots with i ~ j. In this case we 
■ + (1/J) (21) 

QiUj/J) 

and calculate the corresponding converging sum giving the terms in the second line. This anomaly was 
noticed in the Bethe ansatz context in [llj although this phenomenon was known since long in the large N 
matrix integrals or similar character expansions [231 124j . It was proven in [llj to happily cancel in the main 
order of 1/J, even for a more complicated nested Bethe ansatz. 

In our case when J — > oo it is obvious from eqs. (|20H9p that the anomaly docs not contribute to the main 
order and the Bethe ansatz equation becomes a singular integral equation J5J [3] 



27T7lfc =2 

x 



Ctc 



dy go(y) 
x-y 



x € Cfc , k = !,•••, K. 



(22) 



Now we introduce the resolvent and the quasi-momcntum 

1 



G(x) = J2 



1 



xJ — u.j 



p(x) = G(x) + — 



(23) 



as well as the standard definition of density |f| 



p(x) 



1 

2ttI 



(p(x — ie) — p(x + ie j) 



and rewrite it in terms of the following Riemann-Hilbcrt problem 

2p'(x) = 2im k , x&C k , n fe G Z, 



(24) 



(25) 



where rf(x) = \ [p(x + iO) + p(x — iO)] denotes the symmetric (real) part of the quasi-momentum on a cut 
C k , k = l,--- ,K. 

To solve it [7] we notice that p{x) is a function on a hyperelliptic Riemann surface with two sheets related 
by K cuts along the real axis, as shown on the fig[TJ It has a known behavior p(x) ~ Pq + (x) at x = 0, and 
p( x ) ~ S / J + 1 / 2 -f [1/x 2 ) at x — > oo on the first sheet. There are no other singularities. However, as we see 
from eq. (|25p p(x) is not a single valued function but its derivative is. This information is enough to fix p'(x) 
as a single- valued function on hyperelliptic Riemann surface f 2 = Yij=i( x ~ x j) 



a k x 



„k— 1 



The single- valuedness imposes § A dp = 0, I = 1, . . . , K — 1, where Ai, I = 1, • 
surrounding a cut C;. 



(26) 

K are the A-cycles each 



8 These two densities p(x) and q(x) are clearly diferent. By definition q(x) is a smooth function, whereas p{x) is a sum of 
(5-functions. However their 1/J expansion is the same at least for the first two orders. It is also important to note that to obtain 
the 1/J expansion of p(x) one should first expand the rhs of eq,l l24l l in powers of 1/J and then take the limit e — » 0. 
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p(x) for 3-cut solution, S=300 
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Figure 2: Density of roots. The dots correspond to numerical 3-cut solution with total number of Bethe roots 
S = 300 and equal fractions on — 1/6, and ni = { — 1,3,1}. They are fixed from the numerical values of 
the roots by the eg. U6\) . Solid line is the density at J = oo computed analytically from the corresponding 
hyper- elliptic curve, x coordinates of the dots are U3+ ^ J J+1 so that the solitary points in the middle of empty 
cuts are artifacts of this definition. 



The BAE eq. ([25|) become the integer i?-period conditions: 



dp = 2nnj rij = 1, 



(27) 



where the cycle Bj starts at x = oo + on the upper sheet, goes through the cut Cj to the lower sheet and 
ends up at its infinity x = oo_. 

From eq.((2|) we have the quasi-pcriodicity condition for the total momentum expressed by cq.©. In 
some applications, such as the above mentioned spectrum of the integrable matrix of anomalous dimensions 
in Af = 4 SYM theory, we select only purely periodic Bethe states ([7]). This information, together with the 



filling fractions ak = Sk/ J, & = !,••■ , K, such that X)fe=i a k = S/ J, or 



dx p(x) 



fixes completely a solution. For the energy we obtain in this limit 

E = -G' (0). 



(28) 



(29) 



A particular 3-cut solution is demonstrated on the fig. 2 and is compared with the numerical solution of 
exact Bethe ansatz equation. We see that already for 300 roots and J = 600 the description in terms of the 
algebraic curve becomes excellent. 



4 Large J limit and 1/J-corrections from Baxter equation 

Eq. © can be also obtained as the condition that the transfer matrix eigenvalue defining a Bethe state (see 
for example [2"5] ) 

T(u) = W(u + i/2)%^ + W{u - i/2)%=^, (30) 
Q(u) Q(u) 

where Q(u) = Uk=i( u ~ u fc); w ( u ) = u " '■ T ( u ) 

is a polynomial of degree J, which is clear from the very 
construction of a Bethe state in the algebraic Bethe ansatz approach [lj. The Bethe equations ((2|) follow 
immediately from cq. (|30p assuming analyticity (polynomiality) of T(u). 
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Introduce the notations: x = u/J, $(x) = j J2k=i ^°s( x ~ x k), V(x) = logx, 2i(x) = T(Jx)/(Jx) J and 
rewrite (1301) as 



2t(x) = cxp J $ I a; + — j - $(x) + V la; 
Defining the quasi- momentum (now at all orders in 1/J) 

p(x) = $' + V'/2 
and expanding the Baxter equation in 1/ J we get 

1 A/O) V"(s) 



57) - vw 



c.c. 



t(x) 



cosp(x) 
1 



1 - 



sinp(x) 



J \ 2 8 
p"(x) V&(x) 



6 



16 



1 /VO) V"(s) 

2J 2 



1 



(31) 



(32) 



(33) 



J 3 



To find the 1/J corrections to the quasi- momentum, we will expand p(x) = po{x) + jpi(x) + -j?P2(x) + 
0(1/ J 3 ), t(x) = to(x) + jti(x) + -jrtsfa) + 0(1/ J 3 ) and plug it into the last equation. We assume that the 
coefficients of expansion to(x),ti(x),t2(x), . . . are the entire functions on the plane x with no cuts, having 
only an essential singularity at x = 0. 

The quasi-periodicity property of the total momentum reads up to 3 first orders as follows 

p ** = -E^+£i5s? + (1/ j4) = 2nk/J (34) 

and in the purely periodic case we select only the states with k = m J, with integer m. 
4.1 Zero order from Baxter equation 

Let us restore from the Baxter equation the zero order result of the previous section. In the zero order 
approximation we get from eq. (|33[) 

cospo(x) = to (a;) (35) 



or 



2W 



(36) 



As is usual in the finite gap method [26], we expect that since to (a;) is an entire functions all the branch cuts 
of eq.(|26p come from the square root in denominator, after the Bethe roots condense to a set C\, ■ ■ ■ , Ck of 
dense supports in the J — > 00 limit. It is easy to sec from the definition eq.(|32jl that poO) = a+ ^ 2 + (l/^ 2 ) 
when x — > 00 and po(x) ~ 2ttj + (x) when x — > 0. Consequently, we reproduced the general solution (|26H29p 
of 0- 



4.2 1/J correction from Baxter equation 

To find the next, 1/J approximation to the density of roots and to the energy we deduce from eq. (I33|) 



pi = (-p /2 + F"/8)cotp - 



ti 



sinpo 



(37) 



where ti(x) is an entire function on the plane as was mentioned before. We know about po{x) that 

Pq = nrij — nipo, Pq = wrij + nipo, sinpj = — sinp^ (38) 
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and thus we have for the real and imaginary parts of po (x) on the cuts 



(V" \ 1 

™Pi = -z-t - h - — -, = -p cotpo/2. (39) 

V » / sinp 

We will solve these equations below and restore the explicit p\ . 

Since t\(x) is a regular function on the cuts ( [44] ) 

rfx = --np'cottnrp = -pf cotpf /2. (40) 

Moreover we see from cq. (|3~4")) that 

Pi(0)=0 (41) 

and for large x pi(x) should decreases as (1/x 2 ). 

We can write from (|4"0"]) the general solution of this Riemann-Hilbcrt problem 

Pi(») = T^jr, i m^l^M dy + K f a J?L (42) 
4mf{x) J c y{y-x) ~l f( x > 

where f 2 (x) = II j=i ( x ~ x j ) ano - the contour encircles all cuts Ck (but no other singularities) . The first term 
in the r.h.s. represents the Cauchy integral restoring the function from its real part on the cuts and having 
a zero at the origin (the value of the quasi-momentum p{x) at x = 0, oo was already fixed for po) whereas 
the second one is purely imaginary on the cuts, with the polynomial in the numerator chosen in such a way 
that it does not spoil the behavior of p(x) at x = 0, oo^ 

Thus for K < 3 the solution is unique. In particular, for K = 1 we restore from here the 1-cut solution of 
[5] . For K > 3 we have to fix K — 2 parameters aj . To do this we can use K additional conditions ensuring 
the right fractions otj of the roots already chosen for po : 

Px{x)dx = Q, l = l,...,K, (43) 

C, 

in fact only K — 2 of them are linear independent (since we have already fixed the total filling fraction by the 
asymptotic properties of eq. (|4"2")) at x = oo: p\(x) = 0(l/x 2 ). Eq. (|4"T|) also restricts some linear combination 
of the conditions (|43[) V Hence we completely fixed all parameters of our i£f-cut solution for the 1/ J correction 
Pi knowing the zero order solution (algebraic curve) for pq. 

It is also useful to rewrite eq. (14"2"| in the following way 

Pl(») = T^T / fWM^rtv) ^ (44) 

Amf(x) J c Q{y){y-x) 

where Q(x) ~^2^-i ajx^ and the contour of integration encircles all the cuts. Again, aj are fixed by ea. (|43[) . 
Equivalence of this formula to cq. (|42p can be seen from the coincidence of their analytical properties: blowing 
up the contour in the contour integrals of any of these representations we obtain l 10 l 



K-2 



Pl (x) = -p' (x) cotp (x-)/2 + 7 f T £ + J2 JT\> ( 45 ) 

f(x) ^ X n - X ^ f(x) 

where the first term comes from the residue at y = x. Note that only this term contributes to the r.h.s of 
cq. (|4H)) thus showing that we satisfied this equation. The second term exactly cancels all the poles of the first 
term, so that p\ is regular everywhere except the cuts. It is a mcromorphic function in the x plane, having 
no cuts. The last term reflects the freedom of adding K — 2 coefficients before fixing them by conditions 
eq.®. 



9 We could also add terms , jj- , ... but they are too singular at the branch points as we shell see in the next section. 
10 In fact one should take into account an infinite number of residues at y = 0, see [8]. A more regular procedure is to express 
cot as a sum and then do the integration. 
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Figure 3: Quasi-momentum near branch point as a function of the scaling variable v for S = 200. The 
poles corresponds to the positions of Bethe roots Ui- Red dashed line - "exact" numerical value, light grey 
- zero order approximation given by Airy function Ai(a 1//3 x) ; grey - first order and black - second order 
approximation. 



4.3 Equations for 1/J 2 corrections from Baxter relation 

Expanding eq. (f3"3")) up to 1/ J 2 we obtain 



where 



1 1 I2 

P2 = -»d x [cot(po)I] - r-3 - . , v (46) 

2 8x a 2sm(p ) 

= Pl + Pk cotpo . (47) 



sin(po) 2 

We introduced here the notations 

- cospp _ cospo *i cos(2p ) + 5 cos Po . 2 

* X ~ * X + f2 - i2 'l28? + 8?" 24sinp P ° + ^ bo) + ^ (48) 



1 sm po 



so that ii and t 2 are single valued functions. 



Note that above the cut I + = mp\. We will find the explicit solution of these equations later, but we will 
need for that some results of the next section where we study the behavior of p(x) near the branch points. 



5 Double scaling solution near the branch point 

As we stated above the branch point singularities come only from the square roots of the denominator of 
eg. pop . We define an exact branch point as a point where £(:r*) = ±1 . If we approach one of the branch 
points x — > x* we can expand 

t(x) ~ ±[1 - a(x -x,)/2- b{x - x*f/2). (49) 

Note that x*,a,b themselves depend on J. We assume that they have a regular expansion in 1/J (justified 
by self-consistency of further calculations and by numerics) and define x* = Xo + x\j J + . . . . We call xq a 
classical branch point and x\j J a branch point displacement. 

Denoting v = (x — x*)J 2 / 3 which will be our double scaling variable v ~ 1, we get from cq. (|30|) up to 
1/J 2 terms 

/ av bv 2 \ nl , nl W(u + i/2) nl W(u-i/2) , . 
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„—mrvJ 



1/3 



,1/3 



In terms of a new function q(v) = e 
after expansion in 1/J the last eq. takes the form 

1 Avq' + q 1 
In fact, this equation can be easily solved in terms of qo 



Q{x*J + vJ 1 / 3 ), where n is such that t(x*) = e I7m , and 



\2 H y ' 4 



q on 



1 



v 2 1 ( v 4 36 -a 2 



4x 2 ji/a j2/3 l 324 



15a 



1 



1 
J 

126 



9o u 



4ax2j 1 /3 60aJ 2 / 3 



(51) 



(52) 



where ao(v) = Ai(a 1 / 3 ?j) (the Airy function). The second solution of the eg. (f5T|) . Bi(a 1 ^ 3 w) has a wrong 
asymptotic as we will see. The sign oc means that the solution is defined up to a constant multiplier but 
this unknown multiplier doesn't affect the quasi-momentum. Now we can express the quasi-momcntum only 
through our scaling function q(v) 



d v q(v,J) 1 

— -r- ITU A- 

J 2 / 3 J o(w,J)J 1 /3^ ^ 2x* 



1 



- 1 



(53) 



The first two terms in the r.h.s., if we substitute q(v) —* qo(v), represent the principal contribution to the 
double scaling limit near the edge, valid up to the corrections of the order 1/J 2 / 3 . We see from the definition 
(|23p that the zeros of q(v) are nothing but the positions m of Bcthc roots. Thus we know these positions 
with a precision 1/J 2 / 3 (see fig. 3). 

The large v asymptotic will be very helpful in fixing some unknown constant in the 1/J 2 corrections 
given in the next section 

/ \ / \ 



p(x*+vJ 2/3 )=7r?i + — 



JV3 



— y'av— 
1 



1 5 

Av 32v 2 y/av 
-v-' „ ' 



J 2 / s 



/ 



1 



1 



x 2 ^/av I6ax 2 v 2 
V i/J i/J 2 



.,(54) 



where the cut corresponds to negative v for a > 0. Introducing the notation y = vJ 2 / 3 and rearranging the 
terms by the powers 1/J we have 



p{x* + y) 



1 

J 2 



ay 
5 



(a 2 + 12% 3 / 2 
24Va 

1 



1 



1 



4y 8xi^/ay 



a 2 -46 
16a 



(55) 



6-4(a 2 + 126) 
768a; 4 (ay) 3 / 2 



32y 2 y/ay 16ay 2 a; 2 

Doing this re-expansion we assume that J -1 <C y <C 1, trying to sew together the double scaling region with 
the 1/J corrections to the thermodynamical limit. This procedure is similar to the one used in higher orders 
of the WKB approximation in the usual one dimensional quantum mechanics (see for example |27j). 

To compare with p$, pi and pi we have to re-expand around xq 

1 



p{x + y) = p{x* +y) + 



J 2Vv J 2 



Xl 

V 



fax\ 



Wxfajy/Wy 8yy/y 



or, introducing notation 



2A 



x\ 



we get 



p(x + y) 



1 

J 2 



ay 
5 



yfa Ax\a 
(a 2 + 126)y 2 + " 
A 2 



24^Ey 
A 



1 A a 2 -4b 
4y 16a 

6 Jay 



(56) 
(57) 

(58) 



2,2y 2 y/mj 2^/Ey 2 \2y^/mj 64(ay) 3 / 2 768y 2 
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Near the left branch point (i.e for a < and y < 0) we have 



p(xo +y) = nn + 
1 

+ T 2 



(a 2 + \2b)y 2 
5 A 



Z2y 2 Jm) 2^/^ay 2 



1 


' 1 A 


+7 


4y yj— y 


A 2 


b 


2y^Wy 64(ay) 3 / 2 



-46 



16a 
768y 2 



(59) 



+ 



Now we can compare it with our results of the previous sections and fix a, b and x%. 

Let us note that similar Airy type oscillations were observed in the papers on random matrices where this 
behavior occurs near an endpoint of a distribution of eigenvalues [12] . This shows the intrinsic similarity of 
these two seemingly different problems. The large J limit for the spin chain is rather similar to the large size 
N limit in random matrices. We will see this analogy even clearer in the section 7 where we will calculate 
the asymptotics of conserved charges using the results of this section. 

5.1 Comparison with po and p\ 

It is instructive to establish the relations between a, 6, A and the parameters of the algebraic curve which 
completely defines, as we know, p = po + jpi + (jt) up to the first two orders. 

For that we use the expansion (|4"9")l defining a, b and find from eg. ([55)) for y > 



Po{xo + y) = + arccos^o — 



12b 



lay 



2Ay/a~ 



-II 



3/2 



(y 5/ % 



(60) 



in agreement with eqs. (|58l59[l . We can fix a and b up to (1/J) corrections from here through the parameters 
of the solution for po given by eq. ([26|) . 

To calculate a and b up to (1/J) and to fix A, we use the expansion eq. (pl9")) with eg- lpT?]) . Note that we 
have the minus sign in front of ^/ay which ensures the positivity of the density on the cut (i.e. for y < 
and a > 0) p(y) ~ ^Ja{—y)/ir. If we had Bi instead of Ai the sign would be plus and the density would be 
negative. 

Now we compare this near-cut behaviour to p\ . Consider the regular part first 



1 



K,COtpo 



1 

4y 



46 



2" u c 4y 16a 
which is in full agreement with eq. ([55|) . From eq. (|42|) we see that 

Pi (xq + y) - i\ (x + y) ~ — 4 



.3/2 



(61) 



(62) 



where A can be written explicitly, again using the parameters of po given by eq. ([26[) . 
For the example of one-cut solution see eq. ([10111 . 

6 General solution for p 2 and E 2 



Now we have enough of information to construct pi in the most general situation of an arbitrary number of 
cuts. 



We start from a formula which immediately follows from eq. ([46[) 



h = — 2 d x . 



v 

cot(po) ( Pi + y cotpo 



1 



(63) 
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where p\ is given by eq. (|4"2")l . The behaviors near zero and at infinity are the following. Since from eq.© 
and eq.([34]) it follows that G(0) - 2^7tG"(0) = 2nk/J + (^) we can conclude that 

Pa(0) = ^'(0). (64) 

For large x we have again 

p 2 (x) = 0(l/x 2 ). (65) 
Repeating the arguments of the previous subsection we have 

= «75> £ ^ (i? + * [c ° t(i * J ) + '% &y {m) 

where the path C is defined as in eq. (|4"2"l) . Again the first term guarantees that pi satisfies the eq. (l63l) . 
We drop out the ppcothpo f° r simplicity. We can do this since together with f(y) it forms a single-valued 
function without cuts and the integral is given by the poles inside of the path of integration. In fact there 
are only poles at each branch point so that the result can be absorbed into the second term in eq. (1661) . 

So far the second term in cq. (|6"6")) was restricted only by the conditions (j6"4")) and (|6"5|) . Of cause this 
docs not explain why we should restrict ourselves by the fifth power of f(x) in denominator. A natural 
explanation comes from the known behaviour near the branch points (|58I59[) from where we can see that 

p r x i + y ) = | 32F757? - + (a5$5? ~ W^vf 7 * ~ ^S) + (v) ' ai > °' v > ^ 

\ ~:i2y^a iy + + {_2yja iy + 64(a^)3/2 + + (t) ' a t < 0, JJ < 

where all 6K constants eij, b t Ai for i — 1, . . . , 2K are known since they can be determined from the near 
branch point behaviour of po and p\ (|58I59[) . a* and 6, follow from po 




^ " V% + ' «i < °: 2/ < 



(68) 



(69) 



In fact eq. (|67j) gives only two nontrivial conditions for each branch point which are the coefficient before the 
half-integer power of y so that we have 4K conditions. The extra K conditions come from zero A-period 
constraints signifying the absence of corrections to the filling fractions ai. 

p 2 (x)dx = 0, l = l,...,K. (70) 



To reduce the number of unknown constants consider a branch point xq. We can see that for small 
y — x — xq (for simplicity we assume that the cut is on the left i.e. > 0) 

x I fjz) ( 1 \_ 3 A 1 ( b_ _Va\ fl\ 

1 - 4mf(x) % z{z - x) \Az* + 7 ~ l6y 2 V^ ~ 2 V^V 2 + V 3/2 V32a 3 / 2 128 J + \y ) ' 

(71) 

Introduce the following integral 

x I f(z) ( . , Po\ 11/ A 2 36 29JE\ fl\ 

h 55 MR % 7{z—x) [ {P1 + P ° COtPo)Pl COtP0 -12)= 32^ + ^ { V7i - 64^75 + W ) + [y ) 

(72) 
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Figure 4: Resolvent far from branch point as a function of x. Red dashed line - "exact" numerical value for 
one cut solution with S = 10, n = 2, m = 1, light grey - zero order approximation, grey - first order given by 
eg. \4 C 2\) and black - second order approximation given by eg. \7!ty . Note that near branch point (xq = 0.02J the 
approximation does not work and instead of it we should use the Airy function of ea A53\) , like in the usual 
WKB near a turning point. 



so that Ii + 1-2 reproduces the right series expansion near the branch points given by eq. ([M|) and cq. (|55|) . 
Moreover, on the cuts + *0) + I%(x — i0) = since the function under integral is single valued. We can 
simply take 

p 2 (x)=I 1 (x) + I 2 (x)+J^%-\, (73) 

3=0 I{X) 

where the remaining K constants are fixed from eq. ([70]) . Using that p2 (0) = G" (0)/24 we can fix one constant 
Co = G ^24^°^ before imposing the condition (|70p . 

This is our final result for the second quantum correction to the quasi-momentum. In the Appendix B 
we will specify this result for the example of the one-cut solution where it can be made much more explicit. 



7 Energy and higher charges 
7.1 Energy 

To find 1/J corrections to the energy we represent the exact formula eq.{T]) as follows 

I G . (0) + _L C «».(0) + (i r ), ,74, 



E 

J 



where G(x) is defined in eq. (|2"3"|) . We still have to expand G(x) = — +p (x) + -jPi(x) + -jzp^x) + 0(1/ J 3 ). 
Finally, we obtain for the energy: 

E = jE (x) + j^E 1 + j^E 2 + , (75) 

where 

E Q = -G o (0), (76) 
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Figure 5: Relative deviation SE(S)/E(S) of analytical computations of the energy E(S) from its "exact" 
value E exact {S) for the one cut distribution found numerically by Mathematica (solid line corresponds to 
5E(S) = 0), for a finite number of roots S and a finite length J for zero order (light gray), first order (gray) 
and second order (black) approximation. Details are summarized in the table 



# 


1 


2 


3 


4 


5 


m, n 


1,2 


2,1 


1,3 


2,2 


1,5 


Eo 
Ei 
E 2 


12tt 2 
-558.4 
1160 


24ir 2 
-1563 
5464 


16tt 2 
-855.3 
1592 


32tt 2 
-2401 
8982 


24tt 2 
-1563 
1504 


S 
J 
E(S) 
E&(S) 


10 
20 
4.66004 
4.670 


40 
20 
8.54515 
8.619 


7 

21 
5.7359 
5.752 


20 

20 
10.7876 
10.912 


5 

25 
7.0232 
7.070 



and Q(x) = J2^=i bk% k is related to the last term in (|4"2"| . For E 2 we have from cq. (|T3"|) the following 
representation 



= 4 3) (o) _ -i xmfi 



24 ~ P ' 2{0) = 4^7(0)/ + d ^ cot P^- ^ + (Pi + Po Po ) P i cot Po 

ci G%(0)/'(0) 4 3) (0) 
/(0) + 24/(0) 24 ' 



(78) 



Note that for 1-cut we should take ci = 0. We can compare our results with numerical calculations, as it is 
done for a few 1-cut solutions in the figO 



7.2 Local charges 

In this and the next subsection we will calculate local and non-local, or global charges Q r in all powers of 
1/J for the large r from the behavior near the relevant branch point. The idea of this calculation is taken 
from the double scaling approach in matrix models. Namely, we can compare it to the calculation of the 
resolvent of eigenvalues in a gaussian unitary matrix ensemble 

r rt N2 M / N \ oo oo 

/ ^exp ^Tr M 2 Tr(, - M)- 1 ^^~ 2 » ^ ,- 2 "- 1 ff (s ,„ ) (79) 
J [Z7T > V z / g= i n=0 

where M is a hermitian matrix of large size N. The coefficients Hig tTl \ actually give the number of specific 
planar graphs: it is given by the number of surfaces of genus g which can be done from a polygon with 2n 
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edges, by the pairwise gluing of these edges. To extract the large n asymptotics of -ff( g , n ) f° r anv 9 one can 
use that in the large N limit the density (which is the imaginary part of the resolvent on the support of 
eigenvalues) is given by the Wigner's semi-circle law, and the near-edge behavior is described by the Airy 
functional asymptotics [HI HE] showing the traces of individual eigenvalues in the continuous semi-circle 
distribution. We will try to extract the similar asymptotics for the distribution of Bcthe roots. The role of 
1 /N expansion will be played by the 1/ J expansion, whether as the order of the 1 / x expansion in the matrix 
model will be now played by the label r of the charge. 



We start from expanding cq.fK 

oo 

Qr= Yl 



_L (_l)m+l G (r+2m-l)( ) 

jr+2m-i (2m + l)!(r - l)!2 2rn 



(80) 



As we shell see, for large r only the m = term contributes. We express the derivative as a contour 
integral around cuts 

G (n)( Q) = _^LI9W dx . (81) 

For large n only a small neighborhood of the closest to zero branch point xo contributes due to the exponential 
suppression by the l/x" +1 factor. Near the branch point xo we have from eq. (f58|) (see also eqs. (|59!53|l ) 



Gfc(x) = S k0 TTUi - 



1 

2x7 



c k (x-x Q )2 2 \ a \2 2+((x-x Q ) 



i 3ft 

1 Q 



I (-l) k+1 c k (x - x)2 3 2 fc |a| 2 i+Uxo-x) 1 ' 



a > 0, xq < 
a < 0, Xq > 



(82) 



where the universal constants c& can be computed from the known asymptotic of Airy function 



Cfe 



Ai'(z) 



Ai(z) 



Ai(z) 



2^z l / i 



- 

k=0 



4z 3/2 



1 



z 3(n+l)/2 



(83) 



in particular cq — — 1, c\ 



(~2 



32 ' C3 



64' ° 4 



2048' ° 5 



1024' ^6 



414125 
65536 ' 07 



59025 
2048 ' 



These coefficients behave asymptotically as Ck ~ (-l) fe /s! at fc — > oo. 
We assume that /c <C n, r and expand (for xq < 0) 



(y + x )~ n y p dy = \xof +1 ~ n (-l) n i /e- nlog(1 "^dy 



•''(I 



/3+l-n 



(-1)' 



yPe ny dy. 



(84) 



For the last integral the path of integration starts at — oo — iO, encircles the origin in the counterclockwise 
direction, and returns to the point — oo + iO. For the first integral the path is finite: it starts at some point 
— yo — iO where < yo < \xq\ and ends at — yo + iO. The dependence on the yo is exponentially suppressed. 
The last integral is nothing but the Hankel's contour integral 



(v + *o)-Vdy = (-i) n M 



yo 



2ni 



i 



(85) 



similary 



so that 



(y + x )- n (-yfdy = -{xof+^n-^ 1 



4 n) (o) 



(-i)' 



1_ ft 3fc_3 1 _ 

Cftjap 2- K T 5|zq|5 



1 _ fc 3fc _ 3 1 _ 3ft 

/_l\fc+l Cfc|a|2- gj£o|f T 

1 j r(M_i) 



2tti 
(£)). 



a > 0, xo < 
a < 0, Xq > 



(86) 



(87) 
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As we can see from here, only the term with m = in eg. (|80[) contributes at large n. The others are 
suppressed as 1 fn and the final result is 



l_fc 3fc 3 , 



Qk,r 



{-iy Ck ^~%Ll^ 2 ~ (1 + (r~^)) , a > 0, xo < 
where we introduced the notation 



1 oo 1 

Q-jnE^TT- (89) 
Note that Qfc jr is similar to H g n of the matrix model. 

7.3 Non-local charges 

Now we will find the coefficients of the 1 jx n expansion for large n and arbitrary k 

OO -. OO j 

in other words 



Jk x , 

k=o n=l 



d k , n = ^- i £x n - 1 G k (x)dx. (91) 

In fact only the cut with minimal rij contributes for large n, or rather its branch point closest to x — oo. 
The contributions of other branch points are exponentially suppressed. It is enough to consider only a small 
neighborhood of the branch point with maximal |xo,j|. Near the branch point Pk(x) will behave as 

n ( \ a ( 1 A.J c k(x - x )^~^\a\^-i + (x - xo) 1 "*, a > 0, x > , . 

Gfc(x) = dfco Trrii - - — +< . , ; (i as i i , / (i 3fc n r> ( 92 ) 

V 2a; o/ ( (-l) fc+1 c fc (a;o -x)2--\ a \2-2 + ( Xq -x) 1 "^, a < 0, x < 

similarly to the previous section we obtain 



l_k 3k__3 n+ 

Cfc a 2 2 n 2 2x 



1 3k 

2 2 



2 ;\l_\]° (1 + (-- 1 / 2 )) , a > 0, xo > 



(-l) H "- 1 Cia "^^" 2 (1 + K 1/2 )) , a < 0, ,x < 



8 Conclusions 

We showed in this paper on the example of sl(2) Hciscnbcrg spin chain, how to find general solutions to 
quantum intcgrablc problems in a specific thermodynamical limit proposed by Sutherland, and how to find 
various finite size corrections to it. We also propose a double scaling analysis of the near edge distribution of 
Bcthc roots giving some interesting results for the asymptotics of high conserved charges for the finite size 
corrections of any order. 

Our methods borrow some ideas from the matrix models: the size of the chain J is somewhat similar 
to the size of a random matrix N, and the first two l/J finite size corrections to the main limit which we 
calculated are of a similar mathematical nature as the \/N corrections in matrix models. The asymptotics 
of high conserved charges found here in all orders of l/J remind very much the double scaling approach in 
matrix models [T2l [29l [30] . To calculate them we used for that the Airy asymptotics on the edge of the Bcthe 
root distribution, similar to the generic edge behavior in the matrix models. However, it is obvious that by 
finetuning the parameters of the chain and of its large J solutions we can reach various multicritical points 
with a different from Airy classes of universality, probably also similar to those of the multicritical matrix 
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models Q2] or a two cut model [311 [32] ■ Probably, some modern methods of analysis of the doubly scaled 
matrix models, like the Riemann-Hilbert method for the asymptotics of orthogonal polynomials [33| should 
work as well for the quantum intcgrablc chains in this thermodynamical limit. Indeed, it is known that the 
polynomials Qs(u) = Yik=i( u — u k) arc orthogonal for different S with a specific S'-independent measure 

Our methods might be useful for finding the asymptotics of numbers of various combinatorial objects on 
regular lattices related to the intcgrablc models, such the tilting patterns, spanning trees, dimers etc., as its 
analogue was useful for finding the asymptotics of various large planar graphs via the matrix models. 

The methods presented here can be easily carried over to the su(2) quantum chain as well, though some 
peculiarities of this model, like complex distributions of roots and the presence of "string" condensates with 
equally distributed roots [3], should be taken into account. Only slight modifications of our results will 
allow to find the 1/J corrections in the nonlocal intcgrablc deformations of the su(2) spin chain described in 
[351 136j . As for more complicated models solved by nested Bethe ansatz and with the thermodynamical limit 
described by non-hypcrclliptic algebraic curves, the 1/J and 1/J 2 correactions are left to be established. 

In the context of intcgrability in the 4D Yang-Mills theory with M = 4 supersymmetries (SYM) , where 
the hamiltonian of the psu(2,2|4) integrable spin chains describes the matrix of anomalous dimensions, our 
methods could be especially useful. The appropriate thermodynamical limit in terms of non-hyperclliptic 
curve, corresponding to very long operators in SYM was constructed in particular sectors in [H [7[ 1201 137] 
and for the full psu(2, 2|4) SYM chain in [11]. Finding the finite size corrections to these solutions might be 
extremely useful for getting some clues for the quantization of the string sigma model on the other side of 
AdS/CFT correspondence. The finite size correction for the simplest su(3) solution of SYM theory found in 
[38] are not available yet, let alone more complicated solutions and sectors. The finite size corrections to this 
solution, if found, could be compared to the string quantum correction computed in [39] . The corrections 
found here for the most general multi-cut solution in sl(2) sector could shed some light on the structure of 
the first quantum corrections to string solitons describing spinning strings, where for the moment only the 
direct, very cumbersome quasi-classical methods of evaluation of the functional integral near the simplest 
classical solutions are known [15]. 

As for the double scaling limit used here for the refined analysis near the edge of a distribution and 
giving the asymptotics of high conserved charges to any order of 1/J, it is very universal and is not based 
on any particular form of solutions. Near the edge we can forget about the rest of the sheets and cuts in any 
algebraic curve, and the asymptotics will be governed in all cases by the Airy functions. In fact, we think 
that this near edge behavior will persist also in the near-classical regime of the string theory on AdS^ x S 5 , 
and even for some of its recently considered deformations. We expect then that the asymptotics of high 
conserved charges at all orders of 1/J will be also dominated by the same Airy type solution as found here. 
It is already a deeply quantum regime (reminding that of the non-critical string theories in < 2 dimensions) 
and it could be extremely useful as a step to the full quantization and solution of the Metsaev-Tseytlin 
supersgtring and hence of the Af = 4 SYM theory. 
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Appendix A, BAE 1/J expansion 



In this Appendix we give the formulas for the expansion of the rhs of BAE (|2|) in powers of 1/J. The density 
g(x) is defined in cq. (fT6|) . We assume k, S — k ~ J (i.e. far from the ends of the cut) 
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sinh (7rg) 
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7rg'"[coth(7rg)]4 ] + ( j5 ) , 



where we introduced the notation [f(g)] n = /(g) — StTo 1 /' l ' > (0)fr f° r the functions regular at zero. For 
singular functions the Taylor series should be substituted by the Laurent series, so that [f(g)] n is zero for 
g = and has first n — 1 zero derivatives at this point. For example [coth(7rg)] 2 = coth(7rg) 2E£. 



Appendix B, Example: 1-cut 

In this Appendix we express corrections to the energy in terms of infinite sums for the simplest case of 
one-cut solution. For this solution the hyperelliptic curve is a sphere. It is two complex planes connected by 
a single cut. The density of the Bethe roots is given by a simple formula [7] 



2ttx 



We can easily find explicit expressions for aj and bi of eq. ([67|) . With the notation M = y/m(m + n) a,: and 
bj become 



at 



-, b 1 = f 12M + n 2 + 3?i 2 - 47iVAf 2 - 24Af 2 



(V4M 2 + ?i 2 - 2M) 2 3(V4Af 2 + n 2 - 2A-/) 4 V 

(96) 

and 

a ^ = 8Mn 4 7 r 3 = 47r 4 n 6 f 12MV4A/ 2 + n 2 - 3?i 2 + 4tiVA/ 2 + 24M 2 ) 

(V4A/ 2 + n 2 + 2A/) 2 3(V4A/ 2 + n 2 + 2M) 4 V / 

(97) 

It may be more convenient for comparison with string theory results |15j to express A defined by eq. (|58p as 
an infinite sum. We have to evaluate the integral in eci. (|4"!?)) and find A from the behavior near a branch point. 
We compute the integral by poles. To that end we use that the solutions to the equation sin(po(^; ± )) = 
are 

xf = I > 0. (98) 

1 2tt V4A/ 2 +n 2 T V^M 2 + I 2 ~ V ' 
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The points x l=Q are the branch points. They are inside the contour of integration and thus do not contribute. 

1 1 VP + AM 2 1 



Using that f{xf)/xf = ±- and 



X l - Xo.l X; - £0,1 

1 1 



P TTXq 1 

VP + AM' 2 1 

P ^ x 0,2 



(99) 



X{ - X ,2 X[ - X ,2 

We can evaluate the integral eq. (|4"2")l for x — > x (we also take x inside the contour to drop irrelevant 
symmetric part of pi) 

1 I f(v)p'(y) cot p(y) 1 ^ VI 2 + 4M2 
dy -> 2 



2tti 7c 2/(2/ - x) 



1=1 



I 



(100) 



we can conclude that 



V/ 2 + 4M 2 



(101) 



-4i 



I v 



s/l 2 + AM 2 
I 



where the sum should be understood in the zeta-function regularization (a natural explanation why this 
rcgularization gives the right result is given in . A more regular way to express the integral as a sum is to 
expand cot into the sum before integration) 



VP + AM 2 _ 



VP + AM 2 



- 1 "5- 



(102) 



1 = 1 1=1 \ 

We can easily reproduce the result of [8] for E\ in terms of a sum from eq. ([77|) 



Ei = -p'i(0) = 4tt 2 J2 WP + AM 2 



(103) 



1=1 



with f -function regularization assumed. 

We can also express our result for the next correction to the energy Ei given by eq. (|78|) as a double sum. 
We will need the following quantity 



±1 



2n(xffk 



E 

.1=1 



IVP + AM 2 - kVk 2 + AM 2 
P - k 2 



- 1 



Vk 2 + AM 2 1 
2k 2 



Evaluating the integrals in cq. ([78|) wc express E2 as a double sum 

E 2 = -(Xi +I2 + Xs 

where 



/ cotpo) = ~ Vl(0) + £ 



47T»/(0) 



2tt^ ( a/4M 2 + n 2 ± 2 



fc 2 + 2A/ 2 
Vk 2 + 47\/ 5 



(104) 



(105) 



p 1 (x±)-4p' 1 (0) 



4ttz/(0) / 4z 5 



47r 4 M 2 (n 2 + 5M 2 ) 



(106) 
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[ x 0,l 


x 0,2 / 


+ -( 



1 la 



x ,2 \ 96 8a 2 



G' '(0)/'(0) G< 3) (0) 4 



24/(0) 



24 



^M 2 (2n 2 + 11M 2 )tt 4 . 
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Note that in our new notations 1 / Xo,i = 47rM ± 27rV 4M 2 + n 2 . Expressions for at , hi and Ai are given in 
eas. (|96l97j) and eq. (TlOTj) . 
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